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[i]  Application  of  recent  geometric  tools  for  Lagrangian 
coherent  structures  (LCS)  shows  that  material  attraction  in 
geostrophic  velocities  derived  from  altimetry  data  imposed 
an  important  constraint  to  the  motion  of  drifters  from  the 
Grand  Lagrangian  Deployment  (GLAD)  in  the  Gulf  of 
Mexico.  This  material  attraction  is  largely  transparent  to  tra¬ 
ditional  Eulerian  analysis.  Attracting  LCS  acted  as  approx¬ 
imate  centerpieces  for  mesoscale  patterns  formed  hy  the 
drifters.  Persistently  attracting  LCS  cores  emerged  1  week 
before  the  development  of  a  filament  resembling  the  “tiger 
tail”  of  the  Deepwater  Horizon  oil  slick,  thereby  anticipat¬ 
ing  its  formation.  Our  results  suggest  that  the  mesoscale 
circulation  plays  a  significant  role  in  shaping  near-surface 
transport  in  the  Gulf  of  Mexico.  Citation:  Olascoaga,  M.  J., 
et  al.  (2013),  Drifter  motion  in  the  Gulf  of  Mexico  constrained 
by  altimetric  Lagrangian  coherent  stmctures,  Geophys.  Res.  Lett., 
40,  6 1 7 1  -6 1 75 ,  doi:  1 0. 1 002/20 1 3GL058624. 

1.  Introduction 

[2]  Mesoscale  material  patterns  are  commonly  observed 
in  satellite  images  of  the  Gulf  of  Mexico  (GoM).  A 
recent  example  is  the  Deepwater  Horizon  oil  slick’s  “tiger 
tail”  described  by  Olascoaga  and  Haller  [2012].  Despite 
considerable  activity  in  recent  years,  these  material  pat¬ 
terns  remain  poorly  understood,  mainly  because  the  time- 
dependent  flows  that  sustain  them  involve  complex  dynam¬ 
ical  interactions  between  processes  at  both  mesoscales 
and  submesoscales.  Also,  as  these  flows  evolve  rapidly 
in  time,  traditional  Eulerian  methods  often  fail  to  reveal 
material  patterns. 
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[3]  Observations  sufficiently  dense  to  permit  extraction 
of  material  patterns  on  multiple  scales  have  been  rare.  To 
fill  this  void,  more  than  300  satellite-tracked  near-surface 
drifting  buoys  were  deployed  in  the  northeastern  GoM  over 
roughly  1  week  as  part  of  the  Grand  Lagrangian  Deployment 
(GLAD)  in  July  2012.  The  positions  of  the  drifters  from  this 
unprecedented  deployment  were  tracked  for  up  to  6  months. 
(Further  details  relating  to  GLAD  are  provided  in  section  S 1 
of  the  online  supporting  information.) 

[4]  Here  we  show  that  material  attraction  in  geostrophic 
velocities  derived  from  satellite  altimetry  data  exerted  an 
important  constraint  to  GLAD  drifter  motion,  thereby  sug¬ 
gesting  that  near-surface  transport  in  the  GoM  is  signifi¬ 
cantly  influenced  by  the  mesoscale  circulation. 

[5]  To  reach  this  conclusion  we  apply  recent  geometric 
tools  for  Lagrangian  coherent  structures  (LCS) — the  con¬ 
cealed  skeleton  of  material  lines  that  shapes  global  transport 
[Peacock  and  Haller,  2013].  These  tools  enable  objective 
(i.e.,  frame  independent)  and  directly  Lagrangian  detec¬ 
tion  of  LCS  [Haller  and  Beron-Vera,  2012;  Farazmand 
and  Haller,  2013],  and  isolation  of  LCS  cores  [Olascoaga 
and  Haller,  2012].  Attracting  LCS  (generalized  unstable 
manifolds)  form  the  centerpieces  of  material  patterns. 
LCS  cores  (generalized  saddle  points)  are  characterized  by 
persistent  material  attraction  and  represent  precursors  of 
material  instabilities. 

[e]  In  the  analysis  we  consider  daily  mean  positions  of 
GLAD  drifters  in  two  groups  (of  90  drifters  each)  deployed 
4  days  apart  at  nearby  locations,  indicated  by  triangles  in 
Figure  l(left).  Focus  is  placed  on  the  20  days  following 
the  deployment  of  the  first  of  these  groups  (22  July  2012). 
The  large  number  of  drifters  and  the  low  dispersion  during 
this  period  allow  one  to  easily  distinguish  the  formation  of 
patterns  in  their  distributions. 

[7]  Ten  days  prior  to  the  first  deployment  date,  a 
long  chlorophyll  plume  extending  southeastward  from  the 
Mississippi  River  mouth  was  evident,  as  revealed  in  the 
satellite  ocean  color  image  shown  in  Figure  l(left).  Plumes 
like  this  occur  regularly,  possibly  strongly  influenced  by 
the  evolving  mesoscale  field  [Morey  et  al,  2005],  which 
is  dominated  by  the  episodic  shedding  of  Loop  Current 
rings  [Forristall  et  al.,  1992].  The  centerpiece  of  the  chloro¬ 
phyll  plume  is  well  represented  by  attracting  LCS  (black 
curves)  extracted  from  altimetry  using  the  methodology 
described  below.  This  confirms  the  role  attributed  to  the 
mesoscale  circulation.  Connections  between  the  altimetric 
LCS  and  sea  surface  temperature  (SST)  fronts  are  more 
tentative  (Figure  1,  right).  This  is,  to  a  large  extent,  a 
result  of  significant  surface  heating  of  the  Gulf  of  Mexico 
during  the  summer  months,  so  that  SST  is  not  a  good 
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Figure  1.  (left)  Chlorophyll  a  concentration  in  the  northern  Gulf  of  Mexico  on  12  July  2012  derived  from  the  MODIS 
(Moderate  Resolution  Imaging  Spectroradiometer)  sensor  aboard  satellite  Aqua  with  deployment  sites  in  the  two  GLAD 
drifter  groups  considered  indicated  (triangles)  and  altimetric  attracting  LCS  computed  as  30  day  backward  squeezelines 
overlaid  (black  curves),  (right)  Sea  surface  temperature  on  28  July  2012  derived  from  Advanced  Very  High  Resolution 
Radiometer  (AVHRR/3)  sensor  aboard  satellite  NOAA-19  with  LCS  overlaid. 


proxy  for  a  passive  tracer  during  this  time  of  the  year. 
Similar  good  correspondence  is  found  between  altimet¬ 
ric  LCS  and  patterns  acquired  by  GLAD  drifters  in  the 
aforementioned  groups. 

2.  Methods 

[s]  Consider  an  incompressible  two-dimensional  velocity 
held  v(x,  t),  where  position  x  ranges  on  some  domain  of  the 
plane  and  time  t  is  defined  on  the  finite  interval  [a,  /S].  An 
often  used  objective  measure  of  material  deformation  is  the 
Cauchy-Green  strain  tensor, 


C;(xo):=DC;(xo)^DC;(xo),  (1) 


where  Ff^ixo)  :=  x{t;xo,  to)  is  the  advection  map  that  asso¬ 
ciates  times  to  and  t  positions  of  fluid  particles,  which  evolve 
according  to 

X  =  v(x,  t).  (2) 

A  local  normal  growth  measure  of  the  unit  normal,  no,  along 
a  material  line  at  time  to  is  given  by  the  normal  repulsion 
rate  [Haller,  2011]: 


1 

i/no  •  Cl^(xo)-'no 


(3) 


Let  0  <  Ai(xo)  <  A2(xo)  =  Ai(xo)  '  and  ^i(xo)  -L  ^2(^0) 
(dependencies  on  to  and  t  are  omitted  for  simplicity)  be 
eigenvalues  and  (normalized)  eigenvectors  of  (1).  It  follows 
that  p^ixo)  =  a/A2(xo)  >  1  [respectively,  p[/xo)  =  yAi(xo)  < 
1]  for  fluid  particle  trajectories  whose  to  positions  form  a 
material  line  along  which  no  =  ^2(^0)  [respectively,  «o  = 
^i(xo)].  Accordingly,  a  material  line  everywhere  tangent  to 
^i(xo)  [respectively,  ^2(xo)]  is  referred  to  as  a  squeezeline 
[respectively,  stretchline]. 

[9]  Detection  of  attracting  LCS  over  the  time  interval 
[a,  P]  is  achieved  at  any  to  within  this  interval  as  follows. 
Attracting  LCS  a.t  to  =  P  are  identified  with  locally  most 
squeezing  squeezelines  at  to  =  P  obtained  from  backward 
computation  out  to  f  =  a  [Haller  and  Beron-Vera,  2012]. 
Attracting  LCS  at  to  =  “  are  identified  with  locally  most 
stretching  stretchlines  at  to  =  ct  obtained  from  forward 
computation  out  to  t  =  p  [Farazmand  and  Haller,  2013]. 


Attracting  LCS  at  any  a  <  to  <  P  are  obtained  from  back¬ 
ward  [respectively,  forward]  advection  of  attracting  LCS  at 
to  =  P  [respectively,  to  =  a]. 

[10]  LCS  cores  are  segments  of  attracting  LCS  that  expe¬ 
rience  persistent,  as  opposed  to  just  overall,  attraction  over 
the  time  interval  [a,P]  [Olascoaga  and  Haller,  2012]. 
Specifically,  at  to  =  P  an  LCS  core  is  given  by  the  to  = 
P  positions  of  fluid  particles  such  that  (3)  monotonically 
increases  from  to  =  P  to  t  =  a  (for  more  details,  cf.,  section 
S2  of  the  online  supporting  information). 

[11]  The  velocity  held  in  (2)  is  assumed  to  be  of  the  form 

Ax,  t)  =  g/(x2)“‘  V-^  r]{x,  t)  +  V(p(x,  t),  (4) 

where  g  is  the  acceleration  of  gravity;  /(X2)  is  the  latitude- 
dependent  Coriolis  parameter;  r](x,  t)  is  the  sum  of  a  mean 
dynamic  topography  and  the  altimetric  SSH  anomaly  dis¬ 
tributed  by  AVISO  (Archiving,  Validation,  and  Interpreta¬ 
tion  of  Satellite  Oceanographic  data);  and  (p(x,t)  is  such 
that  the  normal  projection  of  v(x,  f)  at  the  coastline  van¬ 
ishes  (details  on  the  altimetric  data  set  and  the  compu¬ 
tation  of  (p{x,t)  are  provided  in  section  S3  of  the  online 
supporting  information). 

[12]  All  integrations  involved  in  the  LCS  analy¬ 
sis  described  below  are  carried  out  using  a  stepsize- 
adapting  fourth-order  Runge-Kutta  method.  Interpolations 
are  obtained  using  a  cubic  scheme.  Differentiation  in  (1) 
is  executed  using  finite  differences  on  an  auxiliary  0.1  km 
width  grid  of  four  points  neighboring  each  point  in  a  regular 
0.5  km  width  grid  covering  the  domain  of  interest. 

[13]  The  GLAD  drifters  followed  the  CODE  (Coastal 
Ocean  Dynamics  Experiment)  design  [Davis,  1985].  The 
drifters  were  drogued  at  1  m  depth  and  tracked  using  GPS 
(Global  Positioning  System)  with  an  accuracy  of  about  5  m 
and  positions  recorded  nominally  every  5  min.  The  raw 
drifter  data  were  treated  to  remove  outliers  and  fill  occa¬ 
sional  gaps,  and  also  low-pass  filtered  with  a  1 5  min  cutoff. 
Our  analysis  uses  daily  mean  drifter  positions. 

3.  Results 

[14]  Figure  2  shows  a  sequence  of  snapshots  of  the  evo¬ 
lution  of  GLAD  drifters  in  the  two  groups  of  drifters  noted 
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Figure  2.  Sequence  of  selected  snapshots  of  the  evolution  of  fluid  patches  as  advected  by  the  altimetry-derived  velocity 
held  (4)  (red  and  blue  areas)  and  corresponding  centerpieces  or  attracting  LCS  obtained  as  the  advected  images  of  forward 
stretchlines  through  the  patches  at  the  initial  time  (red  and  blue  curves).  Also  indicated  are  SI  (red  dots),  S2  (blue  dots) 
drifter  positions,  and  instantaneous  streamlines  of  the  altimetry-derived  velocity  fleld  (gray  curves). 


earlier.  Drifters  in  the  group  deployed  on  22  July  2012, 
denoted  S 1 ,  are  shown  as  red  dots,  while  drifters  in  the  group 
deployed  4  days  later,  denoted  S2,  are  shown  as  blue  dots. 
Several  aspects  of  the  motion  of  the  drifters  are  worth  noting. 
First,  S 1  drifters  take  a  few  days  before  they  begin  to  spread 
and  form  a  discernible  pattern.  This  pattern  is  clearly  evi¬ 
dent  only  after  7  August  2012,  when  a  roughly  100  km  long 
filament  is  observed.  Second,  in  marked  contrast,  S2  drifters 
spread  and  form  a  clear  pattern  rather  quickly  after  deploy¬ 
ment.  Indeed,  by  1  August  2012,  S2  drifters  have  organized 
into  a  long  filament  of  more  than  200  km.  This  resembles 
quite  well  the  “tiger  tail”  Angering  instability  of  the  oil  slick 
during  the  Deepwater  Horizon  spill  [Olascoaga  and  Haller, 
2012;  Huntley  et  a/.,  2011].  Finally,  in  addition  to  organizing 
into  filaments,  both  S 1  and  S2  drifters  exhibit  some  degree 
of  clustering.  This  is  evident,  for  instance,  during  the  forma¬ 
tion  of  the  long  filament  by  S2  drifters,  which  do  not  cover 
the  filament  uniformly. 

[is]  Figure  2  also  shows,  in  gray,  instantaneous  stream¬ 
lines  of  the  altimetry-derived  velocity  fleld  (4)  and,  in  red 
and  blue,  snapshots  of  the  evolution  of  two  patches  of  fluid 
as  advected  by  this  fleld.  Each  patch  is  initially  circular,  cen¬ 
tered  at  the  deployment  site  of  each  drifter  group  on  the 
corresponding  deployment  date.  The  color  convention  is  the 
same  as  that  for  the  drifters.  The  curve  through  each  patch  in 
each  snapshot  is  an  attracting  LCS  obtained  as  the  advected 
image  of  the  forward  stretchline  passing  through  the 


center  of  the  patch  on  the  drifter  deployment  date. 
Stretchlines  were  obtained  from  forward  integration  out  to 
10  August  2012  from  the  respective  deployment  date.  As 
expected,  the  LCS  through  each  patch  forms  its  organizing 
backbone,  providing  an  explanation  for  the  shapes  acquired 
by  the  patches  as  time  progresses.  Identifying  this  from 
visual  inspection  of  the  instantaneous  streamlines  of  the 
altimetry-derived  velocity  fleld  is  not  trivial. 

[i6]  The  extent  to  which  altimetric  LCS  constrain  drifter 
motion  can  be  assessed  by  comparing  the  patterns  acquired 
by  the  drifters  with  those  suggested  by  the  extracted  attract¬ 
ing  LCS.  Quick  inspection  reveals  reasonable  overall  agree¬ 
ment  between  drifter  and  altimetric  LCS  evolution.  Closer 
inspection  reveals  that  the  LCS  initialized  in  the  SI  deploy¬ 
ment  (red  curve)  does  not  form  the  centerpiece  of  the 
distributions  acquired  by  the  drifters  over  the  initial  1 5  days 
of  evolution.  This  LCS  suggests  faster  spreading  than  that 
experienced  by  the  drifters.  Only  during  the  last  5  days  of 
evolution  there  is  a  tendency  of  the  drifters  to  organize  in 
a  direction  suggested  by  the  LCS.  In  marked  contrast,  S2 
drifters  develop  patterns  generally  in  quite  good  agreement 
with  those  suggested  by  the  LCS  initialized  in  the  S2  deploy¬ 
ment  (blue  curve).  This  is  most  clearly  evident  over  the  first 
5  days  or  so  of  evolution,  which  is  characterized  by  the  orga¬ 
nization  of  the  drifters  into  a  long  filament.  The  centerpiece 
of  this  filament  is  very  well  approximated  by  the  LCS.  This 
LCS  caimot  be  claimed  to  form  as  clearly  the  centerpiece  of 
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Figure  3.  Altimetric  LCS  cores  found  within  the  circular  domain  delimited  hy  the  dashed  line  (yellow  dots  with  green 
arrows  attached);  SI  (red  dots)  and  S2  (blue  dots)  drifter  positions;  altimetric  LCS  computed  as  backward  squeezelines 
(light  green  curves);  and  instantaneous  saddle  points  of  the  altimetry-derived  velocity  field  (black  crosses).  The  inset  in  the 
lower  right  panel  illustrates  the  effect  of  a  saddle  point  on  a  material  patch  in  steady  flow.  The  saddle  point,  denoted  p,  is 
the  counterpart  of  LCS  core  and  the  stable  [respectively.,  unstable]  manifold,  denoted  W(p)  [respectively,  is  the 

counterpart  of  repelling  [respectively,  attracting]  LCS. 


the  drifter  patterns  beyond  5  days.  However,  the  drifters  do 
follow  to  some  extent  the  northward  folding  of  the  eastern 
end  of  the  LCS.  This  is  not  true  for  the  southward  folding 
of  the  western  end  of  the  LCS,  which  is  not  accompanied  by 
the  motion  of  the  drifters.  But  the  drifters  do  appear  to  orga¬ 
nize  along  the  red  LCS,  initialized  in  the  SI  deployment. 
Finally,  the  tendency  of  both  SI  and  S2  drifters  to  cluster  is 
not  described  by  the  evolution  of  altimetric  LCS. 

[17]  In  view  of  the  rather  poor  spatiotemporal  reso¬ 
lution  of  the  altimetry  data  set  (with  no  more  than  10 
altimetry  tracks  on  average  per  week  covering  the  entire 
Gulf  of  Mexico),  the  agreement  found  between  altimetric 
LCS  and  drifter  motion  is  remarkable.  The  noted  disagree¬ 
ments  may  well  be  a  consequence  of  the  coarseness  of 
this  data  set,  or  the  dynamics  substantively  deviating  from 
quasigeostrophic  at  the  unresolved  scales.  While  shipboard 
wind  measurements  acquired  during  GLAD  suggest  that  the 
Ekman  drift  is  negligible,  inertial  oscillations  are  known 
to  be  important,  and  other  ageostrophic  motions  such  as 
those  due  to  internal  wave  interaction  with  the  bottom 
topography,  and  mixed  layer  and  frontal  instabilities  may 
be  at  work  too.  Accounting  for  all  these  factors  may  help 
explain  the  differences  in  behavior  between  SI  and  S2 
drifters  (the  latter  were  deployed  over  the  Desoto  Canyon 
on  a  salinity  front).  The  additional  consideration  of  buoy¬ 
ancy  and  finite  size  effects  may  be  needed  to  explain  the 


observed  tendency  of  the  drifters  in  each  group  to  cluster 
[Provenzale,  1999]. 

[is]  We  now  proceed  to  show  that  the  organization  of  S2 
drifters  into  a  long  filament  can  be  anticipated  by  emerging 
LCS  cores  in  the  altimetry  record,  illustrating  the  predictive 
skill  of  LCS-core  analysis  (Figure  3).  This  is  done  by  iden¬ 
tifying  LCS  cores  from  attracting  LCS  at  each  assessment 
time  to,  which  are  computed  as  most  squeezing  squeezelines 
obtained  from  backward  integration  over  a  period  of  1  week. 
The  search  for  developing  LCS  cores  is  restricted  to  a  circu¬ 
lar  region  of  radius  200  km  containing  S2  drifters  (LCS  cores 
outside  this  region  lie  too  far  away  from  the  drifters  to  influ¬ 
ence  their  motion).  Particular  attention  is  paid  to  LCS  cores 
that  persist  over  time  as  these  can  more  effectively  accu¬ 
mulate  momentum  and  hence  more  significantly  affect  the 
motion  of  the  drifters.  Each  LCS  core  foimd  is  represented 
by  a  dot,  indicating  the  position  xl  where  the  backward  nor¬ 
mal  repulsion  rate,  is  maximized.  Attached  to  each 

dot  are  opposing  arrows,  indicating  the  plus  and  minus  back¬ 
ward  squeezing  directions  at  Xq,  ±^i(Xq).  The  magnitude  of 
each  arrow  is  scaled  by  averaged  over  points  xo  on 

the  LCS  core  on  each  side  of  Xg.  The  preferred  direction  of 
stretching  past  the  assessment  time  is  suggested  by  that  of 
the  longest  resulting  arrow. 

[19]  Figure  3  shows  two  LCS  cores  that  can  be  traced 
over  time  from  26  July  2012  out  to  1  August  2012.  These 
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LCS  cores  gradually  loose  strength  over  this  period  of  time, 
with  the  northern  LCS  core  doing  so  at  a  somewhat  faster 
rate  than  the  southern  LCS  core.  Both  LCS  cores  suggest 
asymmetric  material  spreading,  with  the  preferred  direction 
of  spreading  being  southward.  Accordingly,  the  northern 
LCS  core  is  seen  to  be  responsible  for  the  initial  southeast¬ 
ward  motion  of  the  drifters  in  group  S2.  This  southeastward 
spread  is  enhanced  by  the  action  of  the  southern  core  when 
the  drifters  approach  its  area  of  influence,  sometime  between 
28  and  30  July  2012.  This  is  accompanied  by  some  north¬ 
westward  spreading  resulting  from  the  action  of  the  southern 
LCS  core.  By  1  August  2012  the  drifters  in  group  S2 
stretch  out,  delineating  a  long  filament.  Thus,  similar  to 
Olascoaga  and  Haller  [2012],  who  related  the  formation  of 
the  Deepwater  Horizon  oil  slick’s  “tiger  tail”  to  the  emer¬ 
gence  of  an  LCS  core  a  week  before  its  development,  the 
organization  of  S2  drifters  into  this  filament  is  inferred  about 
a  week  in  advance  by  the  emergence  of  the  two  LCS  cores 
discussed  above.  We  emphasize  that  this  inference  relies 
exclusively  on  velocity  data  up  to  the  assessment  time:  no 
velocity  forecast  of  any  kind  is  needed.  We  finally  note  that 
instantaneous  saddle  points  provide  no  indication  about  the 
location,  direction,  or  magnitude  of  the  LCS  cores,  making 
the  conclusions  from  LCS  core  analysis  unreachable  from 
traditional  Eulerian  analysis. 

4.  Conclusions 

[20]  We  have  shown  that  material  attraction  in  geostrophic 
velocities  derived  from  altimetry  data  imposed  an  important 
constraint  to  the  motion  of  GLAD  drifters,  which  suggests 
a  significant  role  for  the  mesoscale  circulation  in  shaping 
near-surface  transport  in  the  GoM.  While  largely  trans¬ 
parent  to  traditional  Eulerian  analysis,  this  conclusion  was 
reached  using  recent  geometric  tools  that  enable  objective 
and  directly  Lagrangian  detection  of  LCS  as  well  as  per¬ 
sistently  attracting  LCS  cores.  Lack  of  resolution  of  the 
altimetry  data  set,  ageostrophic  effects,  and  buoyancy  and 


finite  size  effects  are  among  possible  causes  for  disagree¬ 
ment  between  altimetric  LCS  and  drifter  motion. 
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